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EXTENDED  ABSTRACT 


This  Extended  Abstract  covers  the  activities  on  the  «  Stress  Concentrator  and  Rate  Effect  in 
Formation  of  Adiabatic  Shear  Bands  »  carried  out  during  two  years  (Contract  N°  N68171-95-C- 
9071  and  Contract  N°  N68171-97-C-9003).  The  Report  covers  the  period  from  Nov.  7/1996  to  Jan. 
6/1998.  During  this  period  an  investigation  was  conducted  into  the  effect  of  the  adiabatic  stress 
concentrators  on  development  of  the  adiabatic  shear  bands  (ASB)  and  failure  in  the  shear  mode. 
The  effects  of  the  imposed  velocities,  including  impact  up  to  120  m/s,  on  the  ASB’s  formation  and 
failure  have  been  studied  for  VAR  4340  steel,  ~50  HRC.  In  parallel  to  experiments  an  extended 
numerical  studies  have  been  performed,  up  to  impact  velocities  180  m/s. 

Experiments  on  impact  shearing  have  been  performed  on  four  specimen  geometries  with 
different  stress  concentrators,  from  a  more  uniform  stress  field  (standard  geometry)  to  a  sharp 

notch.  Each  geometry  was  tested  within  a  wide  range  of  velocities,  from  2- 10'^  m/s  to  120  m/s, 
increased  by  19  steps.  At  higher  velocities,  from  10  m/s  to  120  m/s,  a  direct  impact  technique  of  a 
projectile  on  the  Modified  Double  Shear  (MDS)  specimen  have  been  applied. 

For  aU  four  geometries  and  for  all  imposed  velocities  the  results  have  been  obtained  in  the  form 
of  the  nominal  shear  stress  Xn  versus  the  nominal  shear  strain  Fn*  The  characteristic  points  of  those 
curves  like  the  instability  strain  Tnc  and  the  final  localization  strain  Tni  at  failure  have  been  found. 
At  velocities  higher  than  -0.02  m/s  a  substantial  drop  of  the  instability  strain  was  observed.  This 
leads  to  a  substantial  drop  of  the  failure  energy  when  the  impact  velocity  is  increased. 
The  adiabatic  stress  concentrators  facilitate  the  failure  in  the  shear  mode. 

In  order  to  assure  a  well  defined  input  into  the  numerical  studies  a  new  form  of  the  constitutive 
relation  has  been  derived.  The  constitutive  relation,  completely  coupled  with  temperature,  is  based 
to  some  extent  on  the  material  science  approach. 

In  parallel,  an  extensive  numerical  study  have  been  completed  on  wave  effects  and  adiabatic 
stress  concentrators.  The  explicit  FE  code  ABAQUS  was  used  in  this  study.  It  has  been  found 
during  the  first  series  of  calculations  that  the  plastic  work  produces  an  increase  of  the  local 
temperature  around  a  stress  concentrator  which  in  turn  accelerates  failure  due  to  thermal  softening. 
Around  the  notch  tips  a  zone  of  large  plastic  deformation  and  high  temperature  are  formed  very 
rapidly.  The  numerical  study  confirmed  existence  of  the  Critical  Impact  Velocity  in  shear  (CIV). 
It  has  been  found  that  the  trapping  of  plastic  shear  waves  by  adiabatic  instability  is  the  main  cause 
of  the  CIV  in  shear.  The  CIV  can  be  precisely  estimated  by  the  FE  technique.  A  good  correlation 
was  found  for  VAR  4340  steel  between  the  CIV  found  by  FE  code  and  experiment. 

During  the  last  year  of  the  Project  an  introduction  of  a  local  failure  criterion  in  the  FE  code  was 
studied.  The  following  two  local  criteria  were  tested  :  a)  the  local  maximum  inelastic  deformation 
in  the  adiabatic  conditions  ;  b)  the  local  maximum  value  of  the  adiabatic  strain  energy  density. 
Both  criteria  permit  for  numerical  analyses  of  the  fragmentation  processes.  A  general  conclusion 
can  be  drawn  that  in  the  case  of  impact  shearing  the  failure  is  more  sensitive  to  the  adiabatic  stress 
concentrators  at  higher  impact  velocities.  Complete  results  of  the  FE  study  are  given  in  Appendix 
N^^l  to  this  Report. 

List  of  Keywords  :  Adiabatic  shear  bands 
High  strength  steels 
VAR  4340  steel 
Shear  test 
Shear  failure 
Dynamic  plasticity 
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1.  Objectives  of  the  Project  and  general  remarks 

Previous  studies  which  were  carried  out  under  support  of  the  ERO-London  have  shown 
that  the  formation  of  the  adiabatic  shear  bands  (ASB)  in  metals  is  sensitive  to  the  velocity 
of  shearing  [1,  2,  3].  This  is  caused  by  variety  of  effects  like  strain  hardening,  strain  rate 
and  temperature  sensitivity  of  the  flow  stress  and  inertia.  Since  the  imposed  velocity  of 
shearing  can  be  changed  from  very  low  to  low-velocity  impact  (-100  m/s)  and  further  to 
the  ballistic  velocities  (-1000  m/s),  expected  reaction  of  a  material  may  substantially 
differ.  Preliminary  numerical  calculations  by  the  finite  difference  technique  of 
the  adiabatic  shear  localization  have  shown  that  for  a  perfect  layer  in  the  steady  state  of 
shear  (no  stress  concentrators)  the  rate  sensitivity  and  inertia  delay  the  process  of 
localization.  Stages  2  and  3,  [4].  On  the  contrary,  experiments  performed  in 
the  configuration  of  impact  shearing  clearly  demonstrated  that  the  shear  localization  is 
accelerated  when  the  imposed  velocity  is  increased.  This  apparent  contradiction  was 
the  main  reason  to  continue  this  study  on  different  impact  velocities  and  stress 
concentrators  in  adiabatic  shearing.  It  has  appeared  that  the  stress  concentrators  play  a  very 
important  role  in  adiabatic  instabilities  and  adiabatic  failure  of  metallic  materials.  In  order 
to  gain  more  information  on  the  effect  of  stress  concentrators  triggered  in  dynamic 
conditions  (impact)  and  rate  effects  in  formation  of  the  ASB’s,  experimental  studies  have 
been  undertaken.  Different  stress  concentrators  have  been  studied  in  shear  for  VAR  4340 
steel,  -50  HRC. 

On  the  other  hand  numerical  calculations  by  FE  explicit  code  ABAQUS  with  more 
realistic  initial  and  boundary  conditions  have  been  performed  [5,  6].  It  was  intended  to 
verify  available  experimental  results  obtained  in  LPMM-Metz  for  VAR  4340  steel  (-50 
HRC).  The  experimental  technique  of  direct  impact  on  the  MDS  specimen  was  used  during 
testing  of  this  steel,  [1,  7,  8].  The  experimental  finding  was  confirmed  by  the  FE 
calculations  that  the  Critical  Impact  Velocity  in  shear  (CIV)  is  a  kind  of  a  new  material 
constant. 

All  numerical  calculations  for  VAR  4340  steel  have  been  performed  with  a  specially 
developed  constitutive  relations  based  in  part  on  material  science. 

In  order  to  understand  more  the  role  of  the  adiabatic  stress  concentrators  in  shear  failure 
under  impact,  four  geometries  of  the  specimen  with  different  stress  concentrators  were 
assumed.  The  most  acute  stress  concentrator  was  produced  by  the  spark  erosion  technique. 

The  numerical  analyses  were  continued  with  the  FE  explicit  code  ABAQUS  along  with 
the  local  failure  criteria. 

This  approach  permitted  to  study  further  the  effect  of  the  adiabatic  stress  concentrators 
activated  at  different  loading  rates  on  failure.  Two  independent  failure  criteria  were 
introduced  in  the  FE  code.  The  mam  objective  to  resolve  the  discrepancy  between 
experimental  results  (existence  of  the  stress  concentrators)  and  previous  numerical  result 
seems  to  be  accomplished. 
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2.  Experimental  technique  and  Modified  Double  Shear  (MDS)  geometry 

This  experimental  technique,  which  has  been  described  in  details  in  [1,  7,  8],  is 
operational  since  1991  in  LPMM-Metz.  Because  the  test  theory  and  details  of  this 
technique  is  given  elsewhere,  [7],  here  only  a  brief  description  is  offered.  Fig.  1  shows 
the  operational  principles.  Because  of  flexibility  of  specimen  attachment  to  the  loading 
device  (a  hydraulic  machine  or  Hopkinson  tube)  different  geometries  of  the  shear  zone  can 
be  applied. 

The  MDS  specimen  with  its  new  standard  geometry,  which  prevents  plastic  deformation 
and  rotation  of  the  supports,  has  been  discussed  in  [7].  The  deformed  layer  of  the  MDS 
specimen  has  2  mm  and  it  initially  assures  the  uniform  deformation  over  the  gage  length  at 
low  nominal  strain  rates.  The  MDS  specimens  are  loaded  by  direct  impact  using  a  bar 
projectile,  as  it  is  shown  in  Fig.  1.  The  scheme  of  loading  by  direct  impact  permits  for 
greater  flexibility  in  programming  of  the  nominal  strain  rate  in  shear,  simply  by  changing 
the  impact  velocity.  The  rise-time  with  this  type  of  loading  is  very  short  in  comparison,  for 
example,  to  SHTB  technique.  The  flat-ended  projectiles  of  different  lengths  and  diameter 
dp=10  mm,  are  launched  from  an  air  gun  with  predetermined  velocities  V,,  1  m/s  <  Vj  < 
200  m/s.  The  impact  velocity  Vi  is  measured  by  three  sources  of  light  1,  fiber  optics  3, 
three  independent  photodiodes  2  and  two  time  counters  4.  Axial  displacement  6x(t)  of 
the  central  part  of  the  MDS  specimen  is  measured  as  a  function  of  time  by  an  optical,  non- 
contact,  displacement  transducer  8.  The  optical  transducer  reacts  to  the  axial  movements  of 
a  small  black  and  white  target  cemented  to  the  central  part  of  the  MDS  specimen, 
channel  1.  Since  the  impact  end  of  the  projectile  is  black,  the  displacement  of  the  contact 
between  the  projectile  and  MDS  specimen  is  measured  in  the  same  way,  channel  2. 

The  axial  force  transmitted  by  the  specimen  symmetric  supports  can  be  determined  as  a 
function  of  time  from  the  transmitted  longitudinal  elastic  waves  &r(t)  measured  by  strain 
gages  7,  DC  supply  units  11  and  amplifier  10.  All  electric  signals,  that  is  voltages  of 
displacement  Vx(t)  and  Vc(t)  and  transmitted  wave  ej(t)  are  recorded  by  a  digital 
osciloscope  12  and  stored  later  in  the  computer  hard  disc  13  for  further  analyses.  After 
analyses  of  the  recorded  signals  and  elimination  of  time,  a  force-displacement  curve  can  be 

constructed  for  the  MDS  specimen  and  TfF)  and  r(r)  finally  determined,  where  x  and  F 

are  respectively  the  shear  stress  and  shear  strain,  t  is  time  and  F  =  dF  /  dt . 

The  experimental  technique  based  on  the  MDS  specimen  and  direct  impact  has 
appeared  to  be  quite  effective  for  materials  testing  at  high  strain  rates  in  shear,  10^  s'^  <  F 
<  10^  s'^.  Typical  oscillogram  of  such  test  is  shown  in  Fig.  2.  The  MDS  specimen  made  of 
VAR  4340  steel,  ~52  HRC  (standard  geometry)  was  loaded  by  projectile  Lp=200  mm  at 
impact  velocity  Vi=72.3  m/s.  The  two  channel  digital  oscilloscope  with  the  sampling  rate 
1  GHz  was  triggered  by  the  signal  from  the  optical  transducer,  channel  1.  It  is  visible  that 
the  velocity  of  the  central  part  of  the  MDS  specimen  is  almost  constant,  no  vibrations  are 
present.  The  theoretical  time  of  contact  is  tc=2Lp/Co  where  Co  is  the  elastic  wave  speed  in 
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bar,  Co^S.O  inm/|i.s,  finally  tc=80  |as,  the  time  interval  much  longer  in  comparison  to 
the  rupture  time  of  this  specimen,  tf=14  |Xs.  The  transmitted  wave  is  recorded  by  channel  2. 
The  delay  to  propagate  the  elastic  longitudinal  wave  in  the  Hopkinson  tube,  from 
the  specimen  to  the  strain  gage  is  t(i=37.2  [is.  A  special  computer  program  has  been 

developed  to  analyse  the  oscillograms  and  to  determine  T(r)  and  r(r)  curves. 


3.  Effect  of  impact  velocity  on  MDS  specimen,  the  Critical  Impact  Velocity  in  shear 

It  has  been  shown  recently,  [9]  and  [10],  that  during  shear  deformation  imposed  by 
a  high-velocity  impact  plastic  waves  exited  in  a  deformed  material  can  completely  change 
the  mechanics  of  plastic  field.  As  a  rule,  an  intense  plastic  deformation  will  appear  near  the 
impact  end  of  a  specimen.  In  the  case  of  the  MDS  specimen  with  the  gage  length  2.0  mm 
the  nominal  strain  rate  when  the  plastic  waves  start  to  dominate  is  S-IO'^  1/s,  that  is  the 
velocity  of  shearing  is  -100  m/s. 

Since  formulation  of  the  rate  independent  theory  of  elasto-plastic  waves  in  sohds  by 
Karman,  Taylor  and  Rakhmatulin  in  late  forties  and  early  fifties,  for  a  review  see  [1 1],  it  is 
known  that  the  longitudinal  plastic  deformation  can  be  localised  in  thin  bars  by  a  high 
velocity  impact,  [12].  This  deformation  trapping  by  longitudinal  plastic  waves  is  called 
Critical  Impact  Velocity  in  tension,  [12]  and  [13]. 

It  has  been  shown  recently  that  the  Critical  Impact  Velocity  in  shear  can  be  determined 
experimentally  using  the  MDS  specimen  and  the  direct  impact  loading,  [9]  and  [10]. 
An  existence  of  the  CIV  in  shear  was  predicted  by  FE  numerical  method  in  [14].  A  more 
detailed  analytic  study  was  published  elsewhere,  [10],  here  only  a  brief  discussion  is 
offered. 

It  is  clear  that  the  CIV  in  shear  is  closely  related  with  adiabatic  heating  and  thermal 
softening.  The  CIV  in  shear  is  caused  by  an  instantaneous  instability  and  strain  localization 
superimposed  on  plastic  wave  propagation  in  shear.  The  rate-independent  wave 
propagation  theory  has  been  applied  to  analyse  the  CIV  in  shear  and  the  final  formula  was 
derived  in  [10]. 

Preliminary  numerical  calculations  of  the  CIV  for  1018  steel  (French  Standard  XC18) 
have  confirmed  the  usefulness  of  the  analytic  procedure  presented  in  [10].  The  value  of 
CIV  obtained  with  a  simplified  constitutive  relation  was  Vcr=98  m/s  and  rj.=4.9-104  1/s 
for  2.0  mm  gage  length,  the  value  very  close  to  the  one  determined  by  the  MDS  specimen 
and  direct  impact  technique,  Vcr^90  m/s.  The  experimental  technique  based  on  the  MDS 
specimen  and  direct  impact  has  been  applied  so  far  for  the  determination  of  the  CIV  for 
VAR  4340  (52  HRC),  the  value  found  was  =140  m/s,  [10]. 

In  order  to  study  further  the  CIV  phenomenon,  a  numerical  study  of  impact  shearing  of 
a  layer  has  been  performed  by  the  FE  code  ABAQUS.  It  was  intended  to  verify  available 
experimental  results  for  VAR  4340  steel,  -50  HRC,  obtained  by  direct  impact  on  the  MDS 
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specimen,  [1].  In  particular,  the  numerical  study  has  been  focused  on  the  effect  of  the 
impact  velocity  including  the  phenomenon  of  the  Critical  Impact  Velocity  in  shear.  An 
infinite  thin  layer  fixed  at  one  end  and  with  a  small  geometrical  imperfection  in  the  middle 
of  its  height  was  considered.  The  layer  was  submitted  in  the  other  end  to  different  shear 
velocities,  from  quasi-static  to  impact,  up  to  180  m/s.  All  details  of  those  calculations  are 
given  elsewhere,  [5,  6],  the  final  results  are  shown  in  Fig.  3,  [15]. 

Two  modes  of  shear  deformation  of  the  layer  have  been  found:  in  the  first  mode  all 
deformations  concentrates  in  the  middle  of  the  layer,  and  in  the  second  near  the  surface 
where  the  shear  velocity  was  imposed.  The  CIV  could  be  defined  as  a  transition  between 
those  two  modes.  It  has  been  shown  that  impact  velocities  higher  than  the  CIV  lead  to  a 
substantial  reduction  of  the  localization  energy.  Fig.  4. 

All  FE  calculations  have  been  performed  with  constitutive  relations  specially  developed 
to  reflect  behavior  of  VAR  4340  steel.  Thermal  coupling,  strain  hardening  and  rate 
sensitivity  were  accounted  for.  Introduction  of  the  temperature-dependent  strain  hardening 
exponent  has  demonstrated  the  importance  of  the  proper  formulation  of  the  constitutive 
relations  for  all  thermal  instability  problems.  The  constitutive  relations  are  discussed  later 
on  in  this  Report. 

Overall  good  correlation  with  experimental  results  and  physical  intuition  have  been 
achieved. 


4.  Constitutive  relations 


In  the  FE  codes  available  on  the  market  the  constitutive  relations  are  usually  quite 
simple  and  not  specifically  tailored  to  particular  problems.  All  problems  of  plastic 
instabilities  including  thermal  coupling  must  be  analysed  with  a  precise  and  reliable 
constitutive  relation  or  relations.  Precision  of  the  model  reflects  later  the  FE  results.  In 
order  to  approximate  adiabatic  instabilities,  localization  and  wave  effects,  constitutive 
relations  must  include  strain  hardening,  rate  sensitivity  and  temperature  sensitivity  of  the 
flow  stress  within  a  wide  range  of  strains,  strain  rates  and  temperatures.  After  careful 
analyses  of  experimental  data  available  in  open  literature  for  VAR  4340  steel,  ~50  HRC, 
the  following  explicit  form  of  the  constitutive  relation  has  been  worked  out 


x  = 


|I(T) 


1^0 


B 


f  'j’  ^ 


{r„  +  r)” 


+  t. 


^  T  r 

1-— log-^ 

D  *  r 


(1) 


where  B,  |j,o,  t),  n,  m  are  respectively,  the  modulus  of  plasticity,  the  shear  modulus  at 
T=300  K,  the  temperature  index,  the  strain  hardening  exponent  and  the  logarithmic  rate 
sensitivity,  Tq,  Fq  and  are  normalization  constants.  The  temperature  change  of  the  shear 
modulus  is  given  by 
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(2) 


|X(T)  =  Ho[1-AT*-CT*2];  T*  =  T-Trt.  Trt  =  300K 

where  A  and  C  are  constants,  and  T*  =  T  -  is  the  modified  temperature.  In  principle, 
this  version  of  constitutive  equations  can  be  applied  at  RT  and  temperatures  higher  than 
300  K. 

At  the  beginning  a  simplified  version  of  the  constitutive  relation  (1)  was  used  in  the 
preliminary  FE  calculations.  In  this  approach  the  strain  hardening  exponent  n  was  assumed 
as  constant.  It  appeared  that  the  process  of  locahzation  is  unrealistically  extended  into  large 
strains  because  of  a  strong  parabolic  hardening.  It  is  known,  however,  that  the  rate  of  strain 
hardening  in  steel  is  temperature  dependent  and  it  diminishes  with  an  increase  of 
temperature,  [16].  In  the  second  approach  n  was  assumed  as  a  linearly  decreasing  function 
of  the  homologous  temperature,  [16] 

n(T)  =  no(l-TA'J  (3) 

where  Uq  is  the  strain  hardening  exponent  at  T=0  K  and  Tm  is  the  melting  point.  This 
modification  was  found  essential  for  VAR  4340  steel. 

The  structure  of  the  constitutive  relation  has  some  elements  based  on  the  material 
science  approach.  First  of  aU  the  level  of  stress  is  normalised  by  |j,(T)/|j,o  which  takes  into 
account  the  thermal  softening  of  the  crystalline  lattice.  The  first  expression  in  the  brackets 
is  simply  the  internal  stress  and  the  second  one  is  the  rate  and  temperature-dependent 
effective  stress,  [17]. 

The  procedure  of  how  the  constitutive  relations  have  been  developed  is  described  with 
details  in  [5]. 

Total  number  of  constants  in  eqs  (l)-(3)  is  12  and  they  are  given  in  Table  1  and  in 
Appendix  N°1  to  this  Report. 


Table  1.  Material  constants  for  VAR4340  steel,  ~50  HRC 


Const. 

Value 

Unit 

Const. 

Value 

Unit 

B 

1493.0 

MPa 

A 

5.047- 10-4 

1/K 

[lo 

80769.0 

MPa 

C 

1.036-10-3 

1/K2 

0.39 

- 

D 

2914.0 

K 

n 

eq.  3 

- 

T^O 

622.0 

MPa 

no 

0.135 

- 

P 

0.9 

- 

m 

2.776 

- 

X 

38.0 

W/mK 

To 

300.0 

K 

Cl 

7890.0 

kg/m3 

To 

1.6- 10-3 

- 

Cv 

460.0 

J/kgK 

to 

106 

1/s 

Tm 

1812.0 

K 
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Majority  of  results  of  numerical  simulations  of  the  CIV  in  shear  were  obtained  with  the 
modified  version  of  the  constitutive  relation  n(T),  however,  some  results  with  the  constant 
value  of  n  were  also  produced,  mainly,  in  order  to  make  comparison  of  the  FE  results 
where  both  relations  are  used.  Complete  FE  results  are  reported  elsewhere  [5],  in  this 
Report  only  the  most  important  results  are  discussed.  It  should  be  noted  that  modification 
of  the  constitutive  relation  by  n(T)  changes  only  certain  temperature-dependent  terms  in 
eq.  (3),  for  instance,  both  constitutive  relations  have  identical  strain  rate  sensitivity. 

The  most  important  part  of  the  numerical  analyses  is  the  thermal  coupling. 
The  adiabatic  increase  of  temperature  can  be  estimated  from  eqs  (1),  (2)  and  (3).  For  this 
purpose,  the  flow  stress  x  is  taken  as  a  universal  function  of  the  shear  strain  F,  the  shear 
strain  rate  T ,  and  the  temperature  T,  so  that 

t  =  t(r,r,T)  (4) 


It  follows  from  energy  conservation  that 

dT_[i-^(r,T)]T(r,r,T) 
dr~  Po(T)C,(T) 


where  ^(r,T)  is  a  coefficient  representing  stored  energy  in  the  lattice,  po(T)  is  the  density 
of  the  metal,  and  Cv(T)  is  its  specific  heat  at  constant  volume.  The  coefficient  of  stored 
energy  on  average  is  ^=0.05.  Equation  (5)  can  not  be  integrated  in  closed  form  unless 
certain  simplifications  are  made.  Since  po  is  only  weakly  dependent  on  temperature  it  can 
be  taken  as  a  constant;  at  higher  temperature  than  -100  K,  Cy  may  also  be  assumed  as 
a  constant.  Furthermore,  T  can  be  set  equal  to  the  initial  temperature  Tq.  Integration  now 
becomes  possible  so  long  as  one  specifies  r(r) ,  that  is,  so  long  as  one  specifies  changes  in 
strain  rate  during  the  course  of  deformation.  Thus  an  explicit  form  ATCFni,  To)  is  derived 
for  the  increase  in  temperature  per  unit  volume  during  adiabatic  deformation.  With  these 
simplifications,  eq.  (5)  can  be  expressed  as 


AT  = 


1-^ 


p„(TJC,(TJJ 


/t[r,f(r),Tjdr 


(6) 


Once  AT(r,„,To)  has  been  evaluated,  then  eq.  (4)  provides  values  of  At(r,n,To).  Thus 


At  =  T,  (r„„  f,  tJ  -  T  Jr, ,  r,  (t„  +  at)] 


(7) 
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where  tx  and  tA  are  respectively  the  flow  stress  under  isothermal  and  adiabatic  conditions 
of  deformation.  Equation  (7)  can  be  used  to  compare  adiabatic  and  isothermal  stress-strain 
curves  and  to  compare  these  curves  with  those  of  incremental  rate  experiments. 

Thermal  softening  in  shear  due  to  the  adiabatic  heating  leads  to  formation  of  an 
adiabatic  shear  band  (ASB).  The  instability  point  defined  as  (Tm,  Fm)  occurs  at 

(dr A  /  dr)  =  0 ,  Fig.  5.  The  instability  point  can  be  determined  from  equations  (1),  (2)  and 
(3).  It  is  important  to  note  that  the  instability  point  in  adiabatic  shearing  strongly  depends 
on  strain  rate. 

If  the  problem  is  complete,  the  FE  calculations  could  be  carried  out  with  the  thermal 
coupling.  Since  a  large  part  of  the  plastic  work  is  converted  into  heat  the  temperature  of  a 
material  substantially  increases  when  plastic  deformation  advances.  The  balance  of  energy 
with  the  heat  conduction  leads  to  the  following  relation 


n  dT 

dt  ■  dt  ■ 


d^T 

'dy2 


(8) 


where  y  is  the  axis  of  the  heat  conduction,  P  is  the  coefficient  of  energy  conversion,  p=l-^, 
and  X  is  the  heat  conductivity  (Fourier  constant).  When  the  process  is  entirely  adiabatic, 
X=0,  (no  heat  conduction)  the  heating  is  uniform  in  an  elementary  volume.  All  material 
constants  are  given  in  Table  1  and  in  Appendix  N°1  to  this  Report. 

The  thermal  characteristics  which  are  discussed  above  permit  for  numerical  simulations 
of  all  temperature-coupled  problems  as  purely  adiabatic  as  well  as  with  the  heat 
conduction.  The  code  ABAQUS  includes  those  thermal-coupling  problems  in  the  FE 
procedures. 


5.  Adiabatic  stress  concentrators 

Crack  initiation  and  fracture  in  metallic  materials  at  various  loading  rates  have  been  a 
subject  of  research  for  some  time,  for  example  [19].  Since  the  latest  stage  of  ASB  leads  to 
fracture  in  Mode  II,  this  Mode  would  dominate  in  all  cases  of  impact  shearing,  including 
fragmentation  and  perforation.  It  is  assumed  in  the  theory  of  fracture  mechanics  that 
the  crack  tip  has  an  infinitesimal  radius  and  the  stress  concentrator  is  the  most  acute.  If 
the  loading  rate  of  the  small  volume  around  the  crack  tip  is  relatively  short  this  volume  will 
be  heated  adiabatically  and  locally  the  adiabatic  shear  will  be  started.  A  more  detailed 
analyses  of  the  rate  effects  and  stress  fields  at  the  crack  tip  in  Mode  11  is  given  elsewhere 
[20,  21].  The  main  conclusion  after  those  studies  is  that  the  boundaries  of  equal  equivalent 
plastic  strain  at  the  crack  tip  in  Mode  n  can  be  substantially  reduced  when  the  loading  rate 
is  increased.  This  leads  to  a  more  general  conclusion:  the  stress  concentrators  behave  in  a 
more  brittle  manner  at  impact  loading. 
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However,  a  stress  concentrator  may  not  be  acute,  in  such  case  the  plastic  zone  will  be 
larger  and  the  adiabatic  heating  will  be  more  intense.  A  local  failure  criterion  must  be 
applied. 

A  specific  case  is  the  MDS  geometry  loaded  with  a  relatively  high  shear  velocity. 
Propagation  of  the  elastic-plastic  shear  waves  in  such  case  is  non-negligible.  One  example 
with  a  stress  concentrator  (V-type)  is  shown  in  Fig.  6.  In  this  case  the  elastic-plastic  shear 
wave  (adiabatic)  starts  its  propagation  near  the  plane  A-A.  This  wave  is  excited  along  x- 
axis  by  the  shear  velocity  Vq,  and  it  will  propagate  along  y-axis  with  the  wave  speed 
C2(r),  according  to  the  rate-independent  theory  of  wave  propagation,  [11].  The  shear  wave 
shown  as  the  shadow  region  in  Fig.  6  will  approach  the  plane  B-B  with  the  stress 
concentrator,  and  next  the  tip  of  the  stress  concentrator  will  be  deformed  in  the  adiabatic 
conditions  triggering  severe  plastic  deformation.  If  the  carrying  capacity  of  the  material 
will  be  reached  a  failure  will  appear.  It  will  probably  propagate  along  the  maximum  shear 
stress  -  one  of  possible  criteria  of  failure.  This  simplified  picture  may  explain,  though,  an 
interesting  competition  between  two  mechanisms  of  adiabatic  failure.  If  the  impact 
velocity  Vq  is  high  enough  the  shear  wave  can  not  be  propagated  into  the  notch  direction 
and  all  plastic  deformation  will  be  concentrated  near  the  surface  A-A.  This  is  so-called 
Critical  Impact  Velocity  in  shear  (CIV),  [5,  6,  8,  10,  14].  The  CIV  in  shear  is  a  specific 
mode  of  the  adiabatic  failure.  If  the  impact  velocity  is  slightly  lower  than  the  CIV  two 
mechanism  of  failure  will  compete  increasing  the  total  failure  energy  within  a  narrow 
range  of  impact  velocities. 

Both  mechanisms  of  failure  are  sensitive  to  the  adiabatic  stress  concentrators  and 
the  impact  velocity. 


6.  Experiments  with  different  stress  concentrators 

The  main  interest  to  perform  experiments  with  different  stress  concentrators  was  to  study 
the  effect  of  loading  rate  on  formation  of  ASB’s  and  adiabatic  failure.  The  external 
geometry  of  the  MDS  specimens  was  kept  unchanged  and  only  the  shearing  zones  have 
been  modified. 

Four  geometries  were  appUed  as  it  is  shown  in  Fig.  7,  that  is  the  standard  geometry. 
Fig.  7a,  the  U-geometry,  Fig.  7b,  the  V-geometry,  fig.  7c,  and  the  fine  notch  geometry. 
Fig.  7d.  The  ligament  a  as  well  as  the  thickness  b  were  the  same  for  all  geometries, 
a=5.0  mm,  b=6.0  mm  with  exception  of  the  fine  notch  with  a=5.5  mm.  Dimension  of  this 
notch,  cut  by  the  electro-sparking  technique,  was  -0.2  mm.  The  specimens  were  cut  from 
the  leftovers  of  the  cross-rolled  plate  of  VAR  4340  steel  supplied  earlier  by  the  ARO- 
Watertown.  After  thermal  treatment  the  Rockwell  Hardness  was  -50-52  HRC. 

More  details  on  the  experimental  technique,  oscillogram  analyses  and  the  complete  set 
of  results  are  given  in  [15],  the  Report  of  the  end  of  the  first  year. 
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In  this  Report  only  the  mean  values  of  the  following  chai-acteristic  points  are  shown  and 
discussed.  Namely,  the  maximum  shear  stress  tm  and  the  nominal  shear  strain  Fm 
associated  with  Hm,  the  critical  stress  Xc  where  the  fracture  is  detected  after  the  force- 
displacement  plot  and  associated  critical  strain  Fc-  Such  situation  is  shown  in  Fig.  5,  where 
the  point  (Xn,,  Fm)  is  the  instability  point,  the  point  (Xc,  Fc)  is  found  normally  in  Stage  2 
(quasi-brittle  behavior),  or  during  Stage  3  (advanced  adiabatic  shearing).  As  a  rule  Xm  ^  Xc 
and  Fjn  >  Fc.  The  nominal  sheai'  strain  has  been  defined  for  the  same  way  for  ALL  four 
geometries,  Fc=5(t)/h,  where  5(t)  is  the  net  displacement  detennined  ft'orn  the  oscillograms 
as  a  function  of  time  and  h  is  the  gage  length  in  shear',  h=2.0  mm. 

The  fuial  results  of  those  experiments  ai'e  shown  fi'om  Fig.  8  to  Fig.  11  in  the  form  of 

'^m.  Tm  and  Fc  as  a  function  of  the  logarithm  of  the  nominal  strain  rate  F^Cs'^].  In  every 
figure  are  shown  the  mean  values  for  four  geometries  studied.  In  the  range  of  the  small 
strain  rates  all  four  geometries  show  a  positive  rate  sensitivity  of  the  maximum  stress  Xm- 
A  similar'  behavior  is  found  for  the  critical  stress  Xc.  It  appears  that  the  most  “resistant”  is 
the  U-geometry,  obviously  with  the  most  relaxed  stress  concentrator  in  comparison  to  the 
standard  geometi'y.  At  higher  strain  rates  the  negative  effect  of  adiabatic  heating  is 
obvious,  the  stress  Xm  and  Xc  decrease  for  all  four  geometries.  Finally,  at  high  nominal 

sti-ain  rates,  lO^  s-i<F„  <5-104  s'l  (2  m/s<Vi<100  m/s),  values  of  Xm  increase,  and  Xc  is  very 
close  to  Xm  which  indicates  a  quasi-brittle  behavior. 

The  mean  values  of  the  critical  shear  deformations  differ  substantially  for  different 
geometries,  Figs  10  and  11.  Again  the  largest  deformations  are  obsei-ved  for  standai'd  and 
U-geometries. 

Because  of  reduced  plasticity  the  rate  sensitivity  of  shear'  stress  is  also  reduced,  which  is 
visible  in  Fig.  8  and  Fig.  9  for  the  standai'd  geometry.  The  mean  values  of  the  rupture  stress 
is  fi:om  600  MPa  at  low  strain  rates  to  900  MPa  at  strain  rate  -S-lCri  s'L  The  mean  failure 
strain  Fc  increases  within  the  region  of  smaller  strain  rates  from  -0.018  at  r=6-10'4  s'^  up 

to  Fc’=0.03  at  r=1.0  S'^  and  next  diminishes  substantially  to  the  level  -0.01  at  f =10^  s'L 
This  effect  is  due  to  occuixence  of  adiabatic  heating  and  formation  of  ASB’s.  It  must  be 
mentioned,  however,  that  in  overall  the  fracture  behavior  dominates  at  this  level  of 
hardness  (50-«-52  HRC).  A  siiuilai'  behavior  is  found  from  the  U-geometry,  Figs  10  and  11. 

But  the  transition  from  isothermal  to  adiabatic  regimes  at  F  =1.0  s'^  is  more  abrupt. 

Geometries  V  and  with  sharp  notch  behave  as  almost  brittle,  and  the  transition  from 
isothermal  to  adiabatic  regimes  are  bai'ely  visible.  The  mean  value  of  the  maximum  and 
failure  stress  is  about  800  MPa  for  both  cases.  Reduced  volumes  of  plastic  zones  around 
the  notch  tips  reduce  the  nominal  fracture  deformation  determined  at  low  strain  rates  by 
LVTDs.  In  conclusion,  at  this  level  of  hardness  the  level  of  brittleness  obscures  the  typical 
behavior  caused  by  thermally-coupled  plasticity.  It  is  well  known  that  stress  concentrators 
supress  the  volume  of  plastic  zones. 
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The  effect  of  stress  concentrators  and  transition  from  isothermal  to  adiabatic  regimes  of 
failure  is  more  distinct  when  the  plastic  energy  to  failure  is  analysed  as  a  function  of  the 
logarithm  of  strain  rate.  Fig.  1 2  shows  the  plastic  fracture  energies  respectively  for  the  test 
series  from  standard  to  sharp  notch  geometty.  Within  the  domain  of  small  strain  rates 
the  energy  increases  at  increasing  strain  rates,  this  is  the  case  of  standard  and  U- 
georaetries.  The  maximum  level  is  about  18  J  to  19  J  for  the  standard  and  U- geometries. 
The  transition  to  the  adiabatic  regime  is  more  abrupt  for  the  U-geometry,  and  at  high  strain 
rates  the  plastic  energy  to  fractui'e  is  almost  zero.  Such  behavior  is  caused  rather  by  stress 
concenti'ators  and  a  local  adiabatic  instability,  and  not  by,  for  example,  the  CIV  m  shear', 
usually  observed  for  a  more  ductile  metal.  The  stress  concentrators  trigger  fracture-lrke 
processes  at  loading  rates  lower  than  the  CIV  in  shear.  Behavior  of  V-  and  shar'p  notch 
geometries  is  similar.  At  low  strain  rates  the  mean  plastic  energy  to  fracture  is  ar'ound  5  J 
for  both  geometries.  At  high  strain  rates  the  energies  drop  again  to  a  very  low  level. 

A  general  conclusion  emerging  after  all  series  of  tests  on  VAR  4340  steel  (-50- 
52  HRC)  with  different  stress  concentrators  is  that  at  high  loading  rates  this  steel  is  much 
more  prone  to  failure  via  adiabatic  local  failure  than  in  the  quasi-static  conditions. 
Relatively  high  hardness  reduces  substantially  plasticity  and  behavior  is  closer  to  fracture 
mechanics  with  a  reduced  process  zone,  for  example  [22]. 

The  most  important  and  practical  conclusions  on  the  role  of  the  adiabatic  sttess 
concentrators  in  dynamic  failure  can  be  drawn  fi'om  Fig.  12,  where  the  mean  failui'e 
energies  are  shown  as  a  fiinction  of  the  logai'ithm  of  the  nominal  strain  rate  r„.  The  final 
results  resemble  closely  the  ductile  to  brittle  transition  observed  via  firacture  mechanics 
tests  or  Chaipy  test.  When  an  isolated  crack  or  stress  concentrator  is  loaded  at  increasing 
rates  the  total  energy  to  fracture  diminishes,  [19].  The  total  energy  to  failure  differs 
substantially  for  different  stress  concentrators.  The  highest  resistance  is  found,  as  expected, 

for  U-geometry  and  up  to  the  medium  strain  rate,  Fn-l  next  an  abrupt  fall  is  observed 

from  the  level  -20  J  to  -2  J  at  Fn^SO  s"^  a  drop  of  ten  times.  Behavior  of  the  standai'd 
geometry.  Fig.  7a,  is  similar,  the  sharper  stress  concentrators  cause  that  the  total  energy  of 
failure  drops  from  -13  J  at  r„=30  s'^  to  -5  J  at  Fg^lOO  s'F 

The  ti'ansition  occurs  -50  s‘l<fn<-110  s'*,  well  below  the  CIV  value  which  is  ai'ound 
120  m/s,  Fig.  4,  [1,  5].  Behavior  of  two  geometries  with  the  highest  sti'ess  concentrators  is 
similar'.  A  drop  of  energy  is  still  observed,  however  from  the  low  constant  level,  -5  J  at 
-10-3  s-l<f„<-3.0  s-1,  to  -2  J  above  -10  s'K 

In  general,  these  experiments  confii'm  that  the  adiabatic  stress  concentrators  accelerate 
the  failure  processes  failure  at  higher  loading  rates.  In  such  case  the  impact  velocity  may 
be  lower  than  the  CIV  in  shear. 

An  evaluation  of  experimental  findings  by  numerical  methods  is  possible  when  a  simple 
failure  criterion  is  found. 
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7.  Failure  criteria  in  adiabatic  shearing 


Thermal  softening  during  adiabatic  heating  leads  directly  to  instability  and  localization  and 
further  on  to  failure  in  shear.  In  the  case  of  pure  adiabatic  shear  the  condition  for  stability 
is  reduced  to  the  maximum  shear  stress  dx=0.  This  condition  can  be  rewritten  into  the  form 
of  zero  tangent  modules  dx/dr=0  with  the  instability  point  (Xm,  Fm). 

Application  of  the  condition  for  stability  dx/dr=0  and  some  form  of  equation  for 
constitutive  modeling  leads  to  the  final  result  in  the  form  ri„(r,T)  when  the  adiabatic 

process  of  deformation  is  assumed  together  with  a  constant  strain  rate,  [1,  9].  It  has  been 
found  that  an  increase  of  flow  stress  when  the  strain  rate  is  increased  has  a  negative  effect 
on  the  onset  of  adiabatic  instability,  that  is  Tm  is  reduced  when  strain  rate  is  increased. 
The  positive  rate  sensitivity  increases  production  of  plastic  work  converted  into  heat,  this 
process  accelerates  formation  of  the  instability,  so  values  Fm  are  a  diminishing  function  of 
the  nominal  strain  rate.  Several  industrial  steels  tested  so  far  in  shear  at  strain  rates  from 
10'3  to  -10^  1/s  with  the  MDS  method,  [7],  showed  a  substantial  evolution  of  the  stress- 
strain  curve.  Schematic  changes  of  the  x(F)  curves  observed  for  many  industrial  steels  at 
increasing  rates  are  shown  in  Fig.  13. 

In  this  schematic  figure  the  mean  proportions  are  conserved  to  show  the  effect  of  the 
very  high  rate  sensitivity  on  the  yield  stress  and  instability  stress  above  -lO^  1/s.  On  the 
contrary  to  the  positive  rate  sensitivity  of  the  critical  stress  Xm,  the  critical  strains  of 
instability  Fm  always  show  a  tendency  to  diminish  when  the  strain  rate  is  above  -lO^  1/s. 
The  same  happens  with  the  failure  strains.  At  very  high  strain  rates  such  failure  is  even 
more  abrupt,  this  occurs  because  of  a  superposition  of  the  ASB  development  and  plastic 
waves  in  shear. 

In  order  to  introduce  a  failure  criterion  into  a  numerical  code  it  has  been  assumed  that 
the  failure  starts  at  the  instability  point  3x/9F=0  in  the  adiabatic  conditions  of  deformation 
(high  strain  rates).  Thus,  it  leads  to  the  failure  criterion  based  on  the  critical  inelastic 
deformation  Fpc=Fni.  Of  course,  values  of  Fpc  depend  on  strain  rate  and  at  a  low  strain  rate 
(isothermal  conditions)  Fpc— >«>.  At  high  strain  rates  values  of  Fpc  are  substantially  reduced 
according  to  the  scheme  shown  in  Fig.  13. 

Since  all  material  parameters  in  constitutive  relations  (1),  (2)  and  (3)  have  been 
identified  for  VAR  4340  steel,  values  of  Fpc=Fin  were  found.  The  final  result  is  shown  in 

Fig.  14,  where  the  calculated  values  of  Fpc  are  shown  as  a  function  of  logFj,.  As  a 
consequence  of  behavior  shown  in  Fig.  13,  values  of  Fpc  decrease  substantially  up  to  strain 
rate  -10'*  1/s  and  then  they  stabilize. 
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In  general  case  the  failure  rp(.(f)  criterion  leads  to  formulation  in  the  form  of 
the  equivalent  plastic  strain  ip  (ip )  and  equivalent  plastic  strain  rate  ip ,  thus 


V2 


3 


(fill -£22)  +(£22“S33)  +(^33~£ii)  +  ^(^12 +^23 +^31) 


1I/2 


(9) 


with  8p  =  d£p  /  dt  (10) 

when  Ei2?*0  and  £215^0  and  the  remaining  plastic  strain  components  are  zero,  eqs  (9)  and 
(10)  reduce  to 

_  2  j.  2 

ep=;^ei2  and  8p=^£i2  (11) 


and 


(12) 


The  failure  criterion  rp(;(logr)  or  £p  (logSp)  can  be  easily  introduced  in  any  FE  code. 

In  fact,  the  code  ABAQUS  have  already  implemented  such  possibility. 

Another  failure  criterion  which  can  be  derived  using  the  situation  shown  in  Fig.  13  is 
the  criterion  based  on  the  « strain  energy  density »  in  the  adiabatic  conditions  of 
deformation.  Such  criterion  could  be  used  for  a  large  class  of  dynamic  structural  problems. 
The  simplest  failure  criterion  and  most  widely  used  in  quasi-static  metalworking  processes 
has  been  proposed  in  [24], 

C=Jad£  (13) 

0 

where  a  is  equivalent  stress 


a  = 


V2L 


(011-022)  +(^22  “^33)  +(‘^33~<^ll)  +6(012+023+031) 


1/2 


(14) 


£(.  is  the  critical  equivalent  strain,  o^.  is  the  maximum  equivalent  tensile  stress  and  £  is 
the  total  equivalent  strain;  the  constant  C  is  a  material  dependent  critical  value  of  strain 
energy  per  unit  volume.  This  criterion  indicates  that  failure  occurs  when  the  strain  energy 
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per  unit  volume  reaches  a  critical  value.  The  critical  value  is  usually  assumed  as 
a  fracture  strain  obtained  from  tensile  test  either  of  a  rod  or  a  strip.  Failure  is  then  expected 
to  occur  at  the  position  where  the  equivalent  strain  Eg  is  reduced.  In  general,  the  «  strain 
energy  density  criterion  »  is  similar  to  the  critical  strain  criterion.  An  energy  density  failure 
criterion  was  used  for  dynamic  structural  problems  in  [24]  and  also  discussed  in  [25]. 

Since  the  adiabatic  deformation  dominates  in  the  impact  loading  the  failure  criterion 
based  on  the  «  strain  energy  density  »  which  takes  into  account  the  adiabatic  instability  has 
been  specially  developed  during  course  of  this  study  to  be  applied  in  computer  codes.  It  is 
postulated  that  the  adiabatic  failure  will  occur  when  the  total  work  reaches  the  critical 
value  at  Eg,  where  Eg  is  the  equivalent  instability  strain  in  adiabatic  conditions  of 
deformation.  Again,  constitutive  relations  (l)-(3)  were  used  to  calculate  the  energy  at 
different  strain  rates 


Wg  = 


0 


f  =  const. 


(15) 


where  1(0^  is  the  adiabatic  shear  stress  -  shear  strain  relation.  Fig.  5  and  Fig.  13. 

The  final  result  in  the  form  Wg(logr)  is  shown  is  Fig.  14  as  broken  line.  Of  course,  the 
critical  value  of  the  strain  energy  density  diminishes  as  a  function  of  strain  rate  up  to  shear 
strain  rate  -10^  1/s  and  next  it  stabilizes  or  even  slightly  increases.  It  is  interesting  to  note 
that  at  high  strain  rates  10^  s'l<r  <10^  s'^  value  of  Wc  obtained  for  VAR  4340  steel  may 
be  assumed  constant,  Wc«175  MPa.  Since  both  failure  criteria  give  similar  results, 
the  further  numerical  studies  were  continued  with  the  simplest  criterion  based  on 
the  critical  equivalent  strain  in  adiabatic  conditions  of  deformation. 


8.  Numerical  study  of  stress  concentrators  and  adiabatic  failure  in  impact  shearing 

An  extended  series  of  simulations  have  been  carried  out  for  all  four  geometries  of  shear 
layer  with  different  stress  concentrators.  Fig.  7.  The  calculations  have  been  performed  for 
VAR  4340  steel  using  FE  code  ABAQUS/ExpUcit  with  constitutive  relations  (l)-(3)  and 
failure  criterion  based  on  the  ctitical  values  of  deformation  in  adiabatic  shearing. 
The  Appendix  to  this  Report  contains  complete  results  in  the  form  of  tabels  and  figures. 
Here,  only  the  essential  findings  will  be  discussed. 
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In  order  to  compare  experimental  results  with  numerical  simulations  several  shear 
velocities  were  imposed,  from  10  m/s  to  150  m/s.  Every  simulation  for  each  geometry  and 
each  imposed  velocity  Vi  provided  the  following  output  functions; 

(a)  Failure  patterns  (paths  of  fracture)  at  different  times  and  at  different  imposed 
velocities  Vi; 

(b)  Temperature  fields  in  shear  layers  at  the  instant  of  failure,  including  the  hottest  spots; 

(c)  Nominal  shear  strains  at  failure; 

(d)  Displacement  bj  and  82  of  the  central  part  of  Double  Shear  (DS)  geometry  for  all  DS 
geometries,  5^  at  the  beginning  of  failure  (elimination  of  the  first  element)  and  62  - 

displacement  at  the  complete  failure  (elimination  of  the  last  element); 

(e)  Reaction  force  R  determined  at  the  supports  of  each  DS  geometry  for  each  imposed 
velocity  Vi  -  analogy  to  the  measured  force  in  the  MDS  experiments; 

(f)  Determination  of  the  mean  stress  T  via  the  reaction  force  R  as  a  function  of  time; 
an  analogy  to  experiment; 

(g)  The  mean  shear  stress  1:  as  a  function  of  the  nominal  shear  strain  Tn  for  six  different 
velocities  Vi  -  from  10  m/s  to  150  m/s; 

(h)  Energy  of  plastic  deformation  absorbed  for  each  geometry  as  a  function  of 
the  imposed  velocity  Vi. 

This  extended  program  of  FE  simulations  appeared  to  be  very  helpful  to  understand 
better  the  mechanics  of  the  MDS  test  and  to  demonstrate  from  the  other  hand  application 
of  the  developed  failure  criterion. 

For  example,  the  failure  patterns  in  the  final  stage  of  fracture  for  all  four  geometries 
(deformed  configuration)  are  shown  for  three  imposed  velocities  in  figures  from  Fig.  15  to 
Fig.  18.  The  face  B  is  from  the  support  side. 

From  the  fracture  patterns  shown  for  all  four  geometries  in  non-deformed 
configurations  in  Fig.  5  to  Fig.  8  of  the  Appendix  it  is  possible  to  estimate  the  mean 
propagation  speed  of  the  failure,  that  is  ehmination  of  the  consequtive  elements.  It  is 
interesting  to  note  that  the  mean  crack  speed  increases  when  the  stress  concentrator  is 
sharper,  standard  geometry:  a«370m/s;  U-geometry:  -400  m/s;  V-geometry:  -690  m/s; 
Sharp  notch  geometry;  -1700  m/s. 

An  excellent  demonstration  of  the  transition  of  failure  mechanism  for  the  V-geometry  is 
shown  in  Fig.  19,  from  the  failure  starter  at  the  tip  of  the  notch  at  lower  impact  velocities 
(50  m/s  and  130  m/s)  to  the  situation  typical  for  the  CIV  in  shear  (impact  velocities 
160  m/s  and  190  m/s).  This  is  an  analogy  to  discussion  of  Fig.  6  with  the  adiabatic  stress 
concentrator.  This  is  why  superposition  of  plastic  waves  and  stress  concentrators  give 
different  failure  patterns  and  different  external  reactions  at  different  impact  velocities.  For 


22 


example,  calculated  maximum  reaction  force  Rmax  for  the  DS  supports  are  shown  for  all 
four  geometries  as  a  function  of  log(Vi)  in  Fig.  20  (also  Fig.  15  in  the  Appendix).  As 
expected,  when  the  impact  velocity  exceeds  30  m/s  the  reactions  differ  for  each  geometry. 
For  the  standard  and  U-geometry  the  CIV  in  shear  is  well  visible,  but  for  the  V-geometry 
value  of  the  CIV  is  higher  due  to  perturbation  caused  by  the  V-notch.  The  highest 
resistance  to  failure  at  high  impact  velocities  is  observed  for  the  U-notch. 

The  results  of  calculations  plotted  in  the  form  tCFn)  determined  by  the  analyses  of  the 
R(6)  characteristics,  figures  from  Fig.  16  to  Fig.  19  of  the  Appendix,  show  that 
the  begining  of  failure  seems  to  occur  at  relatively  small  nominal  strains,  for  both  low  and 
high  impact  velocities.  The  end  of  fracturing  at  relatively  low  impact  velocities  occurs 
during  descendant  part  of  TCPn)  curves,  but  at  high  velocities  is  closer  to  the  maximum 
stress  determined  via  R(5). 

Analyses  of  the  plastic  energy  expended  to  the  end  of  failure  are  shown  in  Fig.  21  as  a 
function  of  log(Vi),  also  Fig.  13  in  the  Appendix.  The  numerical  results  approach  relatively 
well  the  experimental  results  shown  in  Fig.  12.  At  lower  imposed  velocities  the  standard 
and  U-geometry  show  a  positive  rate  sensitivity  for  both  experiment  and  numerics  with 
maximum  values  of  the  energies  shghtly  lower  than  found  in  experiment.  However,  in 
experiment  the  total  energy  to  failure  was  determined. 

In  both  cases,  numerics  and  experiment,  a  drop  in  energy  is  observed  at  high  velocities, 
with  the  lowest  level  ~3  J.  In  the  case  of  sharp  stress  concentrators  the  level  of  energy  is 
low  also  at  low  impact  velocities,  for  numerics  ~2  J  to  4  J  and  for  experiment  ~6  J  for 
the  V  and  sharp  notches.  On  the  contrary  to  the  experiment  the  V-notch  shows  an  increase 
of  the  plastic  energy  at  high  imposed  velocities,  due  to  the  transition  mechanism  discussed 
earlier.  Perhaps  some  geometrical  and  experimental  imperfections  eliminated  this 
transition  during  experiments. 


9.  Discussion  and  conclusions 

The  experimental  and  numerical  study  of  adiabatic  stress  concentrators  performed  for 
VAR  4340  steel,  ~50  HRC,  have  demonstrated  a  strong  effect  of  the  imposed  velocity  on 
failure  development  in  every  geometry  studied.  At  higher  impact  velocities  a  substantial 
decrease  of  energy  to  fracture  for  VAR  4340  steel  was  revealed.  At  impact  velocities 
around  120  m/s  existence  of  the  CIV  in  shear  has  been  confirmed,  [1,  5,  9],  both  by 
experiments  and  numerics. 

Relatively  high  increase  of  the  local  temperatures  around  the  notch  concentrators 
promotes,  to  some  extent,  a  local  development  of  the  ASB’s.  The  ASB’s  can  propagate 
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later  with  a  high  speed.  In  fact,  formation  of  ASB  is  a  kind  of  fracture  mechanism.  A  more 
detailed  discussion  of  this  subject  can  be  found  elsewhere,  [15]. 

Introduction  to  FE  code  the  failure  criterion  based  on  the  instability  strain  in  plastic 
shearing  has  enabled  to  calculate  fracture  processes  under  impact.  It  has  been  shown  that 
this  local  failure  criterion  is  sensitive  to  the  local  strain  rate.  The  other  possible  criterion  of 
the  adiabatic  failure  can  be  founded  on  the  «  strain  energy  density  »  concept.  Also  in  this 
case  the  criterion  depends  on  the  local  strain  rate.  Both  criteria  lead  to  similar  results  in 
numerical  calculations.  However,  3D  cases  should  be  analyses  further,  including  the  type 
of  element  used  in  FE  codes.  The  preliminary  numerical  analyses  with  the  local  failure 
criterion  confirmed  experimental  findings,  especially  the  reduction  of  failure  energies  at 
high  imposed  velocities.  The  approach  presented  in  this  Project  can  be  applied  in 
numerical  studies  of  more  complex  fragmentation  processes  in  steels. 
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HGURE  CAPTIONS 


Fig.  1  Configuration  of  experimental  setup  for  impact  shearing  of  the  MDS  specimen. 

Fig.  2  OsciUogram  of  the  test  on  MDS  specimen,  VAR  4340  steel  (standard  geometry), 
loaded  by  projectile  Lp=200  mm  with  impact  velocity  72.3  m/s. 

Fig.  3  Nominal  shear  strains  of  instability  Fnc  and  Fni  as  a  function  of  impact 
velocity  Vi. 

Fig.  4  Fracture  energy  versus  impact  velocity  Vi. 

Fig.  5  Schematic  stress-strain  curve  showing  different  stages  of  fast  (adiabatic) 
shearing. 

Fig.  6  Scheme  of  an  adiabatic  stress  concentrator  charged  by  approaching  elastic-plastic 
shear  wave  (shadowed  bar)  excited  by  shear  velocity  Vq. 

Fig.  7  Specimen  geometries  with  different  stress  concentrators. 

Fig.  8  Strain  rate  spectrum  for  the  mean  values  of  Xm,  VAR  4340  steel,  -50  HRC,  four 
geometries. 

Fig.  9  Strain  rate  spectrum  for  the  mean  values  of  Tc>  VAR  4340  steel,  -50  HRC,  four 
geometries. 

Fig.  10  Strain  rate  spectrum  for  the  mean  values  of  Fm,  VAR  4340  steel,  -50  HRC,  four 
geometries. 

Fig.  11  Strain  rate  spectrum  for  the  mean  values  of  Fc,  VAR  4340  steel,  -50  HRC,  four 
geometries. 

Fig.  12  Strain  rate  spectrum  for  the  mean  energy  of  failure,  VAR  4340  steel,  -50  HRC, 
four  geometries. 

Fig.  13  Schematic  evolution  of  stress-strain  curves  of  industrial  steels  at  strain  rates 
higher  than  5-102  i/g;  (i)  ,5.102  i/s;  (2)  -5-103  i/g;  (3)  ,i.io4  i/s;  (4)  ,5.104  i/s; 

Fig.  14  Local  failure  criteria  applied  in  FE  calculations,  critical  shear  strain 
and  critical  energy  density. 

Fig.  15  Final  stages  of  failure  in  deformed  configuration  for  standard  geometry  loaded  by 
three  incident  velocities  Vi,  plane  B  is  the  support  side. 

Fig.  16  Final  stages  of  failure  in  deformed  configuration  for  U-geometry  loaded  by  three 
incident  velocities  Vi,  plane  B  is  the  support  side. 
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Fig.  17  Final  stages  of  failure  in  deformed  configuration  for  V-geometry  loaded  by  three 
incident  velocities  Vi,  plane  B  is  the  support  side. 

Fig.  18  Final  stages  of  failure  in  deformed  configuration  for  sharp-notch  geometry  loaded 
by  three  incident  velocities  Vj,  plane  B  is  the  support  side. 

Fig.  19  Four  patterns  of  failure  for  V-geometry  loaded  at  different  incident  velocities. 

Fig.  20  Maximum  reaction  force  as  a  function  of  impact  velocity;  four  types  of  specimen 
geometry;  logarithmic  scale  of  imposed  velocity. 

Fig.  21  Dissipated  plastic  energy  versus  imposed  velocity;  four  geometries. 
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